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We discuss the semi-classical kinetic equation in the clean limit, with the presence of Berry 
curvatures and magnetic field B, with the aim of applying to Weyl semi-metals. Special attention 
is given to the conservation laws for the collision integrals. It is found that the magneto-resistance 
second order in B is in general negative, with or without Weyl points, though in the later case it is 
in general much smaller. 
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I. INTRODUCTION 

Recently, there are a lot of attentions on magneto¬ 
transport in Weyl senri-metals^— In this paper, I con¬ 
sider the clean limit kinetic equation for an electronic 
system in the presence of Berry curvatures at low mag¬ 
netic fields, in order to better understand transport in 
these systems. There are previous related efforts. There 
is a paper by Son and Spivab^, extending earlier work of 7 
and concentrating on chiral anomaly. They therefore con¬ 
sider the case where there are two (or more) Fermi pock¬ 
ets (valleys), with each pocket enclosing a separate Weyl 
point. Furthermore, they assumed that there is a sepa¬ 
ration of collision time scales, with the intra-pocket scat¬ 
tering time negligible compared with inter-pocket ones. 
They then deduce the magneto-conductivity in this limit 
with the result depending entirely on inter-pocket scat¬ 
tering times. In Kim et aA (see in particular the Sup¬ 
plementary Material), transport equations were derived 
first using a Green’s function approach, then a kinetic 
equation. The first approach is complicated and the re¬ 
sults not easy to use or interpret. In their derivation 
for the kinetic equation, the Fermi exclusion principle 
was ignored, making the validity of results difficult to 
judge. Subsequently, some of the authors ofi provided 
another derivation of the kinetic equation and transport 
coefficients, concentrating on interactions^ and the weak- 
localization corrections^. A related work is also given 
in=, where diffusion equations were derived to study the 
consequences of chiral anomaly. Here we would apply 
a semi-classical kinetic equation, and the clean limit is 
implicitly assumed. We shall consider general ratios be¬ 
tween the intra and inter-valley scattering times, though 
confine ourselves to relaxation time approximations. Par¬ 
ticular attention is made to the conservation of particle 
number under collisions. We then specialize to some spe¬ 
cific cases where the kinetic equations will be solved ana¬ 
lytically. Our results for the electrical conductivities are 


expressed in terms of quantities at the Fermi level. Our 
effort here is rather closely related to^ (but we shall ig¬ 
nore weak-localization corrections which were included 
there). However, seemingly in this work some very spe¬ 
cial approximations have been made, limiting the appli¬ 
cability of their final results. We shall comment on the 
differences as we go along. 

We shall consider a scalar kinetic equation, that is, we 
limit ourselves to the case where the distribution func¬ 
tions are scalars, not matrices. Our formulation involves 
directly only those states near the Fermi level. Hence 
we assume that the electronic bands are not degener¬ 
ate there. Degeneracy at the Weyl nodes away from the 
Fermi level is allowed (and is taken care of by implicit 
connecting conditions as done in*, see also Appendix). 
For simplicity, we shall confine ourselves to zero temper¬ 
ature. 


II. THE KINETIC EQUATION 


Consider then one of the bands. The equation of mo¬ 
tion is given byil 
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where the electric charge is taken to be —e. This sign 
convention is the same as^, but differ from those ir£. 
The rest of the notations are standard, r, k are the 
position and wavevector of the electron. The band is 
specified by the dispersion e k and curvature Q k - E and 
B are the electric and magnetic fields. The solution to 
these equations are: 
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which is also given in£. Eq © and © can also be directly 
read off from the Poisson brackets given inii. Note that 
• B is dimensionless. Eq © and © were derived 
by expansion, order by order in gradients in r and/or 
k. 1 - Therefore rigorously speaking, the equations would 
need modifications near the Weyl nodes, where f2k di¬ 
verges (e.g the phase space volume (1 + -fit B) in eq 
© and Appendix can become negative). We however 
ignore this complication and assume that the equations 
are valid also there, and derive consequences from them. 
It is also useful to note that fik and B are axial vectors, 
while r, k, E are all ordinary vectors. 

The kinetic equation then reads^ (see also the Ap¬ 
pendix), with n(r, k, t) the occupation number, 

d d d 

—n(r, k, t) + r—n(r, k, t) + k—n(r, k, t) = I co i (5) 


where the left hand side accounts for the drift terms in 
phase space, and the right hand side for the collisions, r 
and k from © and © are to be used in this equation. 

We first discuss some properties of eq ©. It can be 
checked that eq © has the solution n(r, k) = /(ek — fj) in 
the absence of an electric field. In this case, using this n 
on the left hand side, the only term left is proportional to 
x B ■ which is also zero due to the vector product. 
On the other hand, the collision integral must also vanish 
for an equilibrium fermi distribution. 

The total density n(v,t) at a given point in space is 
n(r, k, t). The total number current density J is 
]T)k n(r, k,i)r. Explicitly, we have 
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To use the above equations, we should remember that 
the sum over k should be replaced by the integrate (see 
alsoiid£) according to 

E f d 3 k . € , . 

^ (?) 
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Weyl points are characterized by an integer M referred 
sometimes to as the “charge” of the Weyl point. It is 
defined by 
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with the integration is over a neighborhood containing 
the Weyl point. Note that ^ • flk vanishes everywhere 
except at the Weyl point, since 12k is given by the curl of a 
vector potential. Consider the simplest example of a sin¬ 
gle Weyl node where the effective Hamiltonian near this 
point is of the form u*k -r where r are 2x2 Pauli matrices 
and v* > 0. Then the curvature flk, whose zth compo¬ 


nent evaluated via fl z = i where 


ip is the spinor wavefunction (and similarly for x and y), 
is given by flk = for the positive (negative) en¬ 

ergy band, where k the unit vector along k. In this case, 
M = q=l. (If the effective Hamiltonian is rather k • r, 
then M = ±1). Generally, M must still be quantized, 
though flk does not have to have the isotropic form as 
assumed above. It is also helpful below to note that the 
value of ^ ■ flk for the band directly above and below 
the Weyl points must be opposite to each other. 

In the absence of an external electric field, the current 
can be shown to vanish (in contrast to some earlier spec¬ 
ulations, see the Appendix, and alsoi^dZ). To evaluate 
the electric current to linear order in a uniform static 
electric field, we can then subtract eq © from the same 
expression at E = 0, linearizing in E, and then multi¬ 
ply by the electric charge (—e). We then obtain, for the 
electric current, 
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J^ H involves the equilibrium distribution function /(ek — 
p) and integration over all occupied levels. This is the 
familiar anomalous Hall term from the literature^. We 
shall not analyze this further. It is however interesting to 
note that this term can also be written entirely in terms 
of Fermi surface quantities^. J| involves the deviation of 
the distribution function from the equilibrium value, and 
explicitly involves only quantities near the Fermi surface. 
We shall examine it in more detail below. Rigorously 
speaking we need to sum over all bands, but we have left 


out this sum to keep the notations simple. 


III. RELAXATION TIME APPROXIMATION 

To obtain the conductivity tensor relating J| to the 
uniform static field E, we need to evaluate the distribu¬ 
tion function n(k) by solving the kinetic equation (j5j) in 
the uniform steady state 
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In below, we shall consider two special cases and make appr opriate approximations for the collision integral I rn i. 


A. Single Fermi surface 


It turns out that the Berry curvature Ok has interest¬ 
ing implications which do not seem to have been fully 
emphasized in the literature. Let us first examine this 
case and consider a single Fermi surface. (There need 
not be any Weyl points for this subsection). We shall 
assume a simple relaxation time approximation, with 
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where 


5n(k) = n(k) - /(e k - p) (14) 

Here we have already assumed that the local equilibrium 
distribution corresponds also to the chemical potential p 
at equilibrium. This should be contrasted with subsec- 
tion lHI Bl below. Note that then n{ k) = /(ek — p) + 5n(k) 
We have written down a momentum independent relax¬ 
ation time To, though generalization to momentum de¬ 
pendent scattering rates are possible. In the present case, 
the condition of conservation of particle number in colli¬ 
sion implies that 


X/ n (k )= 0 ( 15 ) 

k 

which can be verified a posteriori. On the left-hand-side 
of eq urn there appears a term of the form f (^ x B) • 


^. [To simplify the notation, we shall occasionally leave 
out the subscript k of ek when no confusion would arise.] 
In the absence of Hk, this term is responsible for the 
classical Hall effect. This term makes the solution of eq 
m rather complicated for general S2k- We shall ignore 
this term, thus limiting ourselves to the case where this 
contribution is small compared with the ones that we are 
keeping. (See the comparison below.) 

It is convenient to define 8v{k ) via 
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\de) (1 + fOk-B) 


(16) 


Due to the factor , we need only to consider 8v on 

the Fermi surface, and we have assumed that it can be 
uniquely labeled by k. Generalization to Fermi surfaces 
with more complicated shapes is obvious but would make 
the notations rather clumsy. In this case, we immediately 
get 
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J| now reads, using eq m and C3 
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We convert our results to integration over the Fermi 
surface. Using 
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where = k 2 /[2n 2 \k-VF |] is the density of states per unit 
solid angle dujj, (here v F = (|^) e=jU is the (k dependent) 
Fermi velocity), we finally get 
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An alternative way of writing (flUl) is 
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which we shall also use sometimes below. Here dS; de¬ 
notes an elemental area on the Fermi surface. 

The conductivity tensor er^ can be read off directly. It 
is simply given by 


If B = 0, eq (12211 reduces to the familiar form. The 
effect of the finite field is simply to replace ^ by v, and 
the additional factor 1 + -L2k ■ B in the denominator. 
The subscript a here reminds us that this is only the 
contribution from J|. The anomalous Hall term is needed 
to complete the results. Below we shall analyze some 
special cases of eq 
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where Vj is the jth component of the vector 
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which is just the ordinary expression. The term first order in B is 
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This term is finite only when the system has broken time-reversal symmetry, as expected from general considerations. 
Explicitly, if time-reversal is obeyed, f2k = —fl_k and ^ are odd under k —> — k but the rest of the terms (excluding 

B) are even, so cr^ 1 ] must vanish. The second order terms in B give 
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From this expression, it is immediately clear that each of the diagonal components er u is positive, giving positive 
magneto-conductances. Note that our formulas are valid for a system with finite Berry curvatures even without 
Weyl nodes, so this cannot be assigned to axial anomalies. Explicitly, consider magnetic field along x. Then 
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Here we have simplified our notations and simply write for the x component of fit etc. The off-diagonal components 
of cr b are symmetric under i -n- j. We have 
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These terms may vanish under certain crystal symme¬ 
tries. 

In the above, we have assumed one Fermi surface. 
Weyl points have never been explicitly invoked, and the 
transport properties are deduced entirely by consider the 
distribution functions near the Fermi level. The above 
can trivially be extended to the case of multiple Fermi 
surfaces provided the chemical potential for the Fermi 
surfaces can be taken to be equal, (see next subsection 
when this condition is relaxed). In this case we only 
have to sum over the contributions from each Fermi sur¬ 
face. Eq (l28l) and (l29l) implies that we can have positive 
magneto-conductances even without Weyl points. In the 
present case, this is purely a consequence of the finite 
Berry curvatures. 

Our results are similar to those given in£ but with some 
detailed differences. The comparison is tedious and we 
would not list them in detail here. We comment only 
on some examples. Our (0^1) is consistent with eq (40) 
and (55) in* (when the extra contributions there (cy¬ 
clotron motion and intervalley scattering) are dropped). 
However, they did not reach our (USD-®. For example, 
their (57), (42), when specialized to our present situa¬ 
tion, differ from our K51) and E51) . Seemingly additional 
implicit assumptions were made in2 in the derivations of 
their subsequent equations. 

Let us estimate the order of magnitude of the magneto¬ 
conductance, assuming that the term first order in B 
vanishes. From the above, we see that the ratio of the 
conductivity second order in B relative to the zeroth or¬ 
der (for any direction j ) can be estimated as, e.g., 


where fl is a typical value of f2k on the Fermi surface. 
Here we have also restored the h which was set to unity 
above. A convenient form of this dimensionless ratio can 
also be written as (VL/l 2 B ) 2 , where Ib = ( Tic/eB ) 1/<2 is the 
magnetic length. Numerically, (it is convenient to note 
that ttTic/ e = 2.07 x 10~ 7 in c.g.s. units) we find that, 
with Act = cr(B) — er(0) being the change in conductivity 
by B (we now leave out the subscripts for simplicity), 
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This magneto-conductance is significant only when Q 
is sufficiently large. To have change in conductivity of 
say 2% in a field of one Testla, H has to be larger than 
3 x 10 4 /Bohr 2 . 

In the above, we have ignored the cyclotron motion of 
the electrons (see discussion below eq (fl5l) l which typi¬ 
cally gives a positive magneto-resistance for current flows 
perpendicular to the magnetic field. This contribution 
can be estimated by 

(34) 

where ui c is the cyclotron frequency. The right-hand-side 
of the above can also be written as (l/kpl 2 B ) 2 where l is 
the mean-free path. For our magneto-conductance to be 
more significant than this magneto-resistance contribu¬ 
tion requires kp fl > l. 


B. Multiple Fermi surfaces 
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In this section, we consider multiple Fermi surfaces. 
This situation may include the one where each Fermi 
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surface contains one (or multiple) Weyl point with fi¬ 
nite net charge. A particular limit of this situation has 
already been considered in£, assuming that the intra- 
Ferrni surface scattering is infinitely fast compared with 
inter-Fermi surfaces scattering. (There are additional im¬ 
plicit assumptions taken which we shall discuss further 
below). Our main aim here is to relax this (and the 
other) assumption(s). We shall also indicate how the re¬ 
sults of Ref£ can be recovered as a special limit. (We note 
however that the formulas below would also be applica¬ 
ble to the case where there are simply two Fermi pock¬ 
ets without Weyl nodes, though the resulting magneto¬ 
conductance is expected to be much smaller: see below). 

For definiteness, let us consider two Fermi pockets, la¬ 
beled by 1 and 2, though the generalization to multiple 
pockets is straight-forward with the expense of rapid in¬ 
crease in the number of indices and the number of self- 
consistent equations for the chemical potentials (see (ICT) 
and (|52D below). Let us first consider pocket 1. For the 
intra-pocket scattering, we can write a form analogous to 
eq m 


(ii) _ 5n{ k) 

^ col — 

V) 


(35) 


The superscript in eq (13511 denotes that we are consider¬ 
ing scattering from pocket 1 to 1. 


<5n(k) = n(k) -/(e-/ii) (36) 


now denotes the deviation from local equilibrium. Eq 
(157>1) demands that the relaxation is towards local equi¬ 
librium for this Fermi pocket towards chemical potential 
Hi, which is to be determined later. The condition that 
particle numbers are conserved under intra-pocket scat¬ 
tering implies 

E fe ( k )= 0 (37) 

k61 

similar to eq ED, though in the previous section we can 
directly take Hi to be the equilibrium chemical potential. 

There are also inter-pocket scattering. To write down 
a possible term for the collision integral, it is important 
here to pay special attention to conservation laws2£, in 
particular for particle number. For pocket 1, let us try2i 

41 ? = (™( k ) - /( e - M2 ))/ t 1 2 , (38) 

where H 2 is the chemical potential of pocket 2 and r^ 2 a 
relaxation time, with a similar term for pocket 2. /xi is to 
be determined by (ED, and a similar equation applies to 
H 2 - The conservation of particle number demands that, 
for any distribution n(k), we must have 

•£ " (k) ~~— +E " (k> ~4r ~- = 0 . (39) 

kei Tfe ke 2 Tfc 

This places restrictions on the scattering times r^ 2 and 
t 21 . One possible choice is that both these quantities are 


independent of k. Eq (l39l) then becomes 
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Linearizing in hi ~ M 2 and using (ED and ED, we see that 
this is obeyed if 


Mb) Mb) 

t 12 r 21 v 41 ; 

where we have defined 

Di(B)= [ ^ Pk ( l + -fi k -B) (42) 

Jfcei 47r c 

( with the integral over the Fermi surface 1) and similarly 
for D 2 . It is convenient to define r A so that 


1 1 \J D1D2 

-!2 ~ T x f)^B) 


(43) 


and also with lo2. Here D\ = / kgl is the density 

of states for the Fermi pocket 1 in zero field. It is feasible 
to take r A as field independent, but then r 12 and r 21 
must depend on B. 

Using 


n(k) = [n(k)-f(e- hi)\ +f(t~ Hi) ( 44 ) 


we get 
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This form makes clear that the sum 4<42 ) is finite 

only if Hi 7 ^ 1 * 2 , and represents the decrease (increase) 
in number of particles for Fermi pocket 1 if Hi > (<)/*2 
via intervalley collisions. By construction, the change 
in total particle numbers of the two Fermi pockets via 
collisions must be zero. 

ED , together with the intra-pocket contribution (1351) . 
gives us in total the collision terms 
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where we have defined 

111 \JD 1 D 2 
h t 0 t x Di(B) 


(47) 


fi is in general B dependent. Defining 5v(k) as in eq 
(1161) (note however that 5n(k) is now the deviation from 
local equilibrium, i.e., eq ED), we have 
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With this, eq (1T21) becomes, after dropping the term involving ^ x B under the assumption as stated in the last 
section and linearizing with respect to the electric field, 
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where we have defined 5 hi ,2 = Mi ,2 — M with h the equilibrium chemical potential. (If we have more bands, then more 
terms proportional to chemical potential differences would appear on the right hand side of (HU)). 

For kel, using (l44l) in eq m and using eq (flbl) . we obtain the contribution to J| from Fermi pocket 1 
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It remains to solve for 8 v(k) and <5 mi, 2- The latter are determined by conservation of particle number. Operating 
fkei Pk on ec L 6U and u s in g (EZJ) (which gives f k£l ^p~ k 5v(k) = 0), we get 


f ke i dJ ^Pk eE-& + £(E.B)n k .& 


8 pi — 6^2 = —t 
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Note that we have a similar equation with 1 -o- 2. These 
two equations are consistent with each other only if the 
numerators sum to zero, that is, alternatively, the inte¬ 
gral of the left hand side of (l49|) over k vanishes when 
also summed over all bands. See the Appendix for fur¬ 
ther discussions. 

On the other hand, using (1441) and a similar equation 
for pocket 2, the requirement that the total particle num¬ 
ber be the same as equilibrium requires, with again the 
help of eq (l37l) , 

Di Shi + D 2 6H‘2 = 0 (52) 

Eq (1511) and (l52l) then determines 5 hi, 2 - Once they are 


determined, 5v(k) can then be obtained from eq (H51) and 
the current from (l50l) . 

Let us first consider the special case where the two 
Fermi pockets are located at ±ko and related by in¬ 
version symmetry. In this case, if k € 1 lies near ko, 
then — k £ 2 lies near — ko, with fl(— k) = O(k) but 
^|_k = — Tjjj| k . Note that these two Fermi pockets to¬ 
gether must enclose zero total charge M. Eq (15)21) then 
implies that 5p 1 = — 8 p 2 , as expected, and so the left 
hand side of (1511) simply becomes 25Hi, with furthermore 
D\ = D 2 (and Di = D 2 ) which we shall simply denoted 
as D (£>). Eq (H51) then gives 


5is(k) = —fi 
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Inserting this into eq (1501) gives us finally, using the short-hand v defined in eq (1231) for the contribution from pocket 
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Pocket 2 gives an equal contribution, so that the total J| is just twice of J| x . The first term of (15H) is the same as 
what had in the last section (eq (l20l) h except that To —► fi, which is in general also field dependent. The second term 
is a new contribution. ,_ 


The zero field conductances and finite field magneto- conductances can be obtained by expanding eq (l54l) . (re- 






























call that fi, D , v are all field dependent). The result¬ 
ing formulas for the general situation are quite lengthy. 
We simply indicate here the relations between our results 
and that in Ref^. Our expression for the current general¬ 
izes that of— in that we have also included the effects of 
deviation from local equilibrium 8v(k ), which gives the 
first term in eq (1541) (and also a correction to the second 
term). This contribution can be dropped if To << r x . 
With this approximation, we now have 
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In this equation, f k£1 in general vanishes, leav¬ 

ing us only a second order term in B. Defining 
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we get, for the total current from both pockets, (with 
Mi = M 2 = M in the present case) 


a transport property. Using © and we have 
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where n = ^/| ^ | is the unit vector normal to the Fermi 
surface. Employing Gauss theorem, Mi is then seen to 
be proportional to the total outward flux of fl(k) from 
the Fermi surface. Since fi(k) is divergence free, Mi is 
thus quantitized. If there is a single Weyl point inside the 
Fermi surface, M is just M, the charge of the Weyl point 
it encloses. More generally, it is a sum over the charges. 
If one assumes that only one isotropic Weyl point of unit 
charge is enclosed, using D = k 2 F / (2 tt 2 vf) = p 2 /2ir 2 v F 
with kp the Fermi momentum, then eq (1571) reduces to 
that of— except with a difference in notations. (t x /2 —> r 
there: see eq (|38l) L 


,e,(2) 

' a,tot 




e 4 M 2 
(47t 2 ) 2 c 2 D 


B(B E) 


(57) 


The magneto-conductance depends on M, an integral on 
the Fermi surface, as one expects on general grounds for 


Let us return to the more general case where the two 
Fermi pockets are not related by symmetry. Solving for 
from eq m and (l52l) and substituting in eq (l50l) . 
we find that the current is given by the same expres¬ 
sion as in (1541) . provided we use the replacement 


J 


1 

2D 


2 


r x Dj 


d 2 


(TLA ,) 1 / 2 d 1 + d 2 


1 - 


Ti (DiD 2 ) 1/2 (^i 

r x b x b 2 



in the second term. In the tq <C t a limit, eq (l57l) holds 
with the replacement D —► ( D\D 2 ) 4 / 2 in the denomina¬ 
tor. 

Despite the more involved form of (l57l) . the order of 
magnitude for the magneto-conductance can still be es¬ 
timated in similar manner as in subsection IIII Al Since 
the zero field conductivity is dictated rather by the first 
term in there is an extra factor ~ t x /tq for the 

ratio Acr/cr(0) in eq (l33l) . If the Fermi pockets do not 
enclose any Weyl points, M = 0 and the second term of 
eq (l54l) vanishes, the estimate of the last section applies 
giving much smaller magneto-conductances. 


IV. SUMMARY 


showed how our results reduce to those^ in a special limit. 
We hope that our more generalized equations can provide 
better understanding of transport in Weyl systems, espe¬ 
cially when system parameters are known in more detail. 
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In this paper, employing a semi-classical kinetic equa¬ 
tion with relaxation time approximation and taking into 
account the effect of Berry curvature in Weyl metals, we 
evaluate the magneto-conductivity at low fields in the 
clean limit. We pay special attention to conservation of 
particles when writing down the collision integrals. We 


Appendix A: Further discussions on the kinetic 
equation 

a. Let us denote the phase space volume by 

7£(r,k, t). In our case, 1Z = 1 + |Hk • B. On general 
grounds, one expects 1Z to obey a continuity equation in 
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phase space 

< 97 ? r) <9 

i* + s-^ + sr (kK > = ° < A1 > 

This equation has also been given ii>^ though in a slightly 
different form. The kinetic equation, in the absence of 
collisions, can also be written as 

d(nlZ) d d • 

“5T + 55' (rnR) + 5k' (k,lE) = 0 <A2) 

The continuity equation (1T21) follows from using (IA1D in 
(IA2I) (after adding the collision contribution). 

Let us check examine dm) directly by using © and 
Q, in the case where the material parameters such as 
fj k are time and space independent. We have 

dn _ e <9B 
dt c dt k 


0 6 dc 

— • (r n) = efi k • (V x E) + -(— ■ n k )V ■ B 

and 

J--(k7e) = -(E-B)(v k -n k ) 

a k c 

For a general point k, since the curvature f2k is given 
by the curl of a vector potential, |;-fl k vanishes. Hence, 
in the absence of Weyl points, eq (IA1I) is indeed valid after 
we invoke the Maxwell equations V • B = 0 and + 
V x E = 0. With a Weyl point, there is an additional 
source term ^ ■ Jlk- 

This has some unusual consequences. The continuity 
equation has to be replaced by£ one with an extra source 
term S , he., 

+ = (A3) 


(disappearing to) the origin and going towards finite k. 
The total number is precisely given by the term S. The 
presence of this term should not alarm us: for the band 
right below, M must be exactly opposite, meaning simply 
that the additional particles from one band must come 
from the other through this Weyl point. 

Alternatively, one can consider the change in the ki¬ 
netic energy of a particle due to the drift terms. This 
is given by k • Besides the ordinary terms propor¬ 
tional to E • ^ describing the work done by the electric 

field, we have the additional term — ^r(E- B)Ok • ^ 
Hence, if E ■ B < 0 say, then the kinetic energy is in¬ 
creasing if fik and the velocity are antiparallel, show¬ 
ing that the energy of the particle would increase with 
time. If the particle is originally below and near the Weyl 
point, it would eventually goes through the Weyl point 
and emerges above it. 

For 1Z, this source term can be understood in the same 
manner as the source term for extra particles discussed 
two paragraphs above. We note however that our deriva¬ 
tion of magneto-conductance in text never invoked ex¬ 
plicitly the continuity equation (IA3I) . However, the dif¬ 
ference between (n times) (IA1I) and (IA2I) yields the term 
proportional to (E • B)llk • f in, e.g., m and eventu¬ 
ally responsible for the appearance of M in eq (1561) and 
the mangeto-conductance (ED. 

b. We next show here that, (i) in the absence 
of an external electric field the electric current vanishes, 
and (ii), the eq ED is consistent with itself with lf>2. 
Equation © gives, with n(k) given by the equilibrium 
Fermi function f(e k — /r), the expression 




d 3 k 




/(ek~M) (A5) 


with a sum over all bands, where these bands need not 
intersect the Fermi level. We note that the right hand 
side of (1511) involves a very similar integral as in eq (IA5I) . 
Let us define 


where 

S = -^ rc (E-B)f(e*-ii)M (A4) 

The right hand side of (IA3I) denotes the contribution to 
the change in density from collision, but the term S on 
the left denotes the contribution entirely from the drift 
in phase space. 

This source term in the continuity equation has the 
same origin as the chiral anomaly discussed in~ for strong 
magnetic field. In our case, it can be understood by 
examining ©• To simplify our discussion, let us pre¬ 
tend that the Weyl point at k = 0 and the dispersion is 
isotropic around this point. We then have = fjv- For 
the band with energy above the Weyl point, £k = e* +v*k 
for some v* > 0. Then we have for k, up to the factor 
1 + §n k • B, the term ^(E ■ B)O k = ' B )#- 

If (E ■ B )M < (>)0, then we have particles coming from 



d 3 k 

(2^)3 


dek e de k 
dk + c dk 


f2 k )B 


S(e 


Eq (1511) is then 


e') 

(A6) 


- 6fi 2 


_x f i(aO-E 
(D 1 D 2 ) 1/2 


(A7) 


and the equilibrium current is 


J = 


de / F tot (e')/(e' - n) 


(AS) 


where F to t is the sum J] Fj over bands j. 

It is then sufficient to prove that the sum Ftot(e / ) van¬ 
ishes for any e'. Let us call the first and second contribu¬ 
tions in eq m F“ and Fj, with similar definitions for 
F^. The vanishing of F“ ot is easiest to see. Using OTT) . 
we see that the x t component of F“ ot is proportional to 
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the the surface integral of the dot product of Xi with the 
Fermi surface normal n ( = ^ /1 ^ I) and hence vanishes 
since the Fermi surfaces must form close surfaces. 

Let us now consider F* ot . The contribution from band 
1 is proportional to f kel with the integral 

now over a surface at energy e'. We show that, when 
summed over all bands, this quantity vanishes. One way 
to proceed is to note that the above integral is propor¬ 
tional to Mi defined in (1551) 1. and then uses the Fermion 
doubling theorem^. However, we may also proceed as 
follows. We write this sum, using the inverse of eq (TT9l) . 
as — f ( 2^)3 Ok • fjL where the integral is now over 
all momentum k, and the sum is over all bands, and 
n = /(e — e'). Using integration by parts and periodic¬ 
ity of n(k)J7k over the Brillouin zone, this sum can be 
rewritten as 



d 3 k d n 
J2^ n dk ' k ' 


It is obviously zero when all bands are empty. Its value 
can change only when a point where ^ ■ L2k ^ 0 changes 
its occupation, but then the contribution above and be¬ 
low this Weyl point cancels, since we can replace n by 
unity and the bands above and below the Weyl points 
have opposite ^ ■ flk- [It is may be helpful also to 
put the argument slightly differently: (a) the sum of 
fk over the bands vanish when e ’ is below 

all Weyl points using the integration by parts argument 
given above; (b) This sum does not change with e’ when 
the Weyl points are not crossed (also via the arguments 
above); and finally (c) when e 1 11 increases from slightly be¬ 
low to slightly above a Weyl point, such that say band 1 
changes from a small hole pocket to an electron pocket, 
the integral f kl ^ for this band can be veri¬ 

fied to be unchanged by considering a Hamiltonian h ■ r 
with h linear in k.] This completes the proof, provided 
all bands are accounted for. 
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